Continued fraction representation of the Coulomb Green's operator 
and unified description of bound, resonant and scattering states 
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If a quantum mechanical Hamiltonian has an infinite symmetric tridiagonal (Jacobi) matrix form 
in some discrete Hilbert-space basis representation, then its Green's operator can be constructed in 
terms of a continued fraction. As an illustrative example we discuss the Coulomb Green's operator 
in Coulomb-Sturmian basis representation. Based on this representation, a quantum mechanical 
approximation method for solving Lippmann-Schwinger integral equations can be established, which 
is equally applicable for bound-, resonant- and scattering-state problems with free and Coulombic 
asymptotics as well. The performance of this technique is illustrated with a detailed investigation 
of a nuclear potential describing the interaction of two a particles. 



PACS number(s): 03.65.Ge, 02.30.Rz, 02.30. Lt 
I. INTRODUCTION 

Green's operators play a central role in theoretical 
physics, especially in quantum mechanics, as they ap- 
pear in fundamental equations governing the dynamics of 
physical systems. The formalisms based on Green's op- 
erators represent alternative, but essentially equivalent 
way of describing the same systems as methods based 
on differential equations. However, since Green's opera- 
tors are related to the integral equation formalism, their 
use is more advantageous in many cases than that of 
traditional differential approaches: boundary conditions, 
for example, are automatically incorporated in the for- 
malism. Although the concept of Green's operators fre- 
quently appears in basic textbooks on the level of fun- 
damental equations, in practical calculations their use 
is usually avoided, and various approximations to the 
Schrodinger equation are preferred instead. The rea- 
son certainly is that the evaluation of Green's operators 
is much more complicated than the direct treatment of 
the Hamiltonian using standard tools of theoretical and 
mathematical physics. 

The discrete Hilbert-space basis representation of the 
Green's operator is very advantageous since it makes 
the solution of the Lippmann-Schwinger-type integral 
equations possible on an easy-to-apply, yet very general 
way. If the integral equation possesses good mathemat- 
ical properties, then the potential term can be well ap- 
proximated on a finite subset of the Hilbert-space basis 
and the integral equation can be solved without any fur- 
ther approximation. This concept has already been elab- 
orated in several different forms. In Refs. harmonic 
oscillator functions were used, which allowed the repre- 
sentation of the free Green's operator, and thus the cor- 
responding approximation method can handle problems 
with free asymptotics. In Refs. || || Coulomb-Sturmian 



basis was applied, which allowed the inclusion of the long- 
range Coulomb potential into the Green's operator, and 
thus led to the exact treatment of the Coulomb asymp- 
totics. This latter approach has also been extended to 
solving the three-body Coulomb problem in the Fad- 
deev approach. So far good results have been reached 
for bound-state and below-breakup scattering-state 
problems || . These results have showed the efficiency of 
the discrete Hilbert-space expansion method in solving 
fundamental integral equations. We note that in all the 
previous approaches [|]-|| the Coulomb Green's matrix 
has been determined in a rather complicated way with 
the aid of special functions. 

In Ref. H we have presented a rather general and 
easy-to-apply method for calculating the discrete Hilbert- 
space basis representation of the Green's operators of 
those Hamiltonians, which have infinite symmetric tridi- 
agonal (i.e. Jacobi) matrix forms. The procedure ne- 
cessitates the evaluation of the Hamiltonian matrix on 
this basis, a rather common task in quantum mechan- 
ics. Then the elements of the Jacobi matrix are used as 
input in the calculation of the Green's matrix in terms 
of a continued fraction. This way of calculating Green's 
matrices simplifies the calculations considerably and the 
representation via continued fraction provides a readily 
computable way. The combination of this new way of cal- 
culating the Coulomb Green's matrix with the technique 
of solving integral equations in discrete Hilbert-space- 
basis representation results a quantum mechanical ap- 
proximation method which is rather general in the sense 
that it is equally applicable to solving bound-, resonant- 
and scattering-state problems with practically any po- 
tential of physical relevance. And all this is provided at 
a very little cost: in practice only matrix elements of the 
Hamiltonian are required. 

The structure of this paper is the following. In sec- 
tion |n] we consider the Coulomb Green's operator in 
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Coulomb-Sturmian representation. We show that this 
Coulomb Green's matrix can be given in terms of a con- 
tinued fraction. In section [II the solution method of 



with coefficients 



the Lippmann-Schwinger equation in Coulomb-Sturmian 
space representation is explained, and the method is ap- 
plied to describe bound, resonance and scattering solu- 
tions of a model nuclear potential. Finally conclusions 
are drawn in section IV. 



II. CONTINUED FRACTION REPRESENTATION 
OF THE COULOMB GREEN'S OPERATOR 

Let us consider the radial Coulomb Hamiltonian 



H 



2m V dr 2 



1(1 + 1) 



where m, I, e and Z stand for the mass, angular mo- 
mentum, electron charge and charge number, respec- 
tively. The Coulomb-Sturmian (CS) functions, the 
Sturm-Liouville solutions of the Hamiltonian (2.1) [Q, 
appear as 



(r\n) 



(n + 2/ + l)! 



exp( 



-br)(2br) l+1 L {2l+1 \2br) 



(2.2) 



where & is a scale parameter, n is the radial quan- 
tum number and denotes the generalized Laguerre 
polynomials |l3f| . With the biorthonormal partner de- 
fined by (r|n) = {r\n}/r these functions form a dis- 
crete basis. We denote the Coulomb Green's operator 
as Gf(z) = (z — Hp)' 1 and we consider here its CS 
matrix elements G^ n , = (n\Gf\n'). 

The starting point in this procedure is the observa- 
tion the the matrix J_^ n , — (n\(z — Hp)\n') possesses an 
infinite symmetric tridiagonal i.e. Jacobi structure, 



J^ n = 2(n + l + l)(k 2 -b 2 



Arab 



- Z^ 



and 



[n(n + 2l + l)] 1/2 (fc 2 + b 2 ) 
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(2.3) 



(2.4) 



where k — (2mz/h 2 ) 1 ^ 2 is the wave number. The main 
result of Ref. || is that for Jacobi matrix systems the 
N'th leading submatrix G^j N ^ of the infinite Green's 
matrix can be determined by the elements of the Jacobi 
matrix 



n C(N) 
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where C is a continued fraction 

UN 



C= — 
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(2.5) 



(2.6) 
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* = -l?,i/J?,i + i • (2-7) 



In Rcf. [|| we have shown that the continued fraction 
C, as it stands, is convergent only for negative energies, 
but it can be continued analytically to the whole complex 
energy plane. Since the Ui and di coefficients possess the 
limit properties 



u — lim Ui = 

i — >oc 
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and 



d= lim d, = 2(k 2 -b 2 )/(k 2 + b 2 ) , 



the continued fraction appears as 

u N 
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(2.8) 



(2.9) 



(2.10) 
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The tail w of C satisfies the implicit relation 

u 



d + w ' 

which is solved by 

w± = {b±ik) 2 /(b 2 + k 2 ) 



(2.11) 



(2.12) 



Replacing the tail of the continued fraction by its explicit 
analytical form w±, we can speed up the convergence 
and, which is more important, turn a non-convergent con- 
tinued fraction into a convergent one JlT[ . In fact, by us- 
ing w± instead of the non-converging tail we perform an 
analytic continuation. In Ref. Q we have shown that w+ 
provides an analytic continuation of the Green's matrix 
to the physical, while w_ to the unphysical Ricmann- 
sheet. This way Eq. (2.6) together with (2.5) provide the 
CS basis representation of the Coulomb Green's opera- 
tor on the whole complex energy plane. We note here 
t hat with the choice of Z — the Coulomb Hamiltonian 
( |2.l| ) reduces to the kinetic energy operator and our for- 
mulas provide the CS basis representation of the Green's 
operator of the free particle as well. 

We can immediately test G c by calculating eigenval- 
ues belonging to attractive Coulomb interactions. Fig. 
I] shows \G (E)\. The poles coincide with the exact 
Coulomb energy levels up to machine accuracy. We stress 
that, from the point of view of determining the energy 
eigenvalues, the rank N of the matrix and the specific 
choice of the CS basis parameter are irrelevant. Any 
representation of the Coulomb Green's operator exhibits 
all the properties of the system and our Green's matrix 
contains all the infinitely many eigenvalues. This is espe- 
cially interesting if we compare with the usual procedure 
of calculating eigenvalues of a finite Hamiltonian matrix 
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which could only provide an upper limit for the N low- 
est eigenvalues. Our procedure does not truncate the 
Coulomb Hamiltonian, because all the higher J_ nn , ma- 
trix elements are implicitly contained in the continued 
fraction. We note that G has already been calculated 
before |3-|| by using %Fi hypergeometric functions jL3[ . 
In Ref. B] we have shown that the two formalism lead to 
numerically identical results for all energies and the con- 
tinued fraction representation possesses all the analytic 
properties of G c also in practice. 



but they turned out to be also very efficient in solving 
integral equations in discrete Hilbert space basis repre- 
sentation |l6| . The choice of 

N l-cxp{-[a(n-jV-l)/(jV + l)] 2 } 
a * = l-exp(-«2) (3J) 

with a ~ 5 has proved to be appropriate in practical 
calculations. 

Let us introduce an approximation of the potential op- 
erator 



JV 



III. SOLUTION OF THE 
LIPPMANN-SCHWINGER INTEGRAL 
EQUATION IN CS REPRESENTATION 



In this section, after Refs. [|3|-|6|, we recapitulate the 
solution of the Lippmann-Schwinger equation in CS rep- 
resentation. We suppose that the Hamiltonian Hi is split 
into two terms 



Hi 



Vi 



(3.1) 



where Vi is the asymptotically irrelevant short-range 
potential. The Green's operator of Hi is defined by 
Gi(z) = (z — Hi) -1 , and it is connected to Gf(z) via 
the resolvent relation 



G l {z) = G?{z) + GY(z)V l Gi{z) 



(3.2) 



The wave function l 1 ?!^) describing a scattering process 
satisfies the inhomogeneous Lippmann-Schwinger equa- 
tion 



\M ±) ) = M ±) )+Gf(E±iO)V l \^ ±) ) , 



(3.3) 



where |3>; ) is the solution of the Hamiltonian Hp 
possessing scattering asymptotics. The bound- and 
resonant-state wave functions satisfy the homogeneous 
Lippmann-Schwinger equation 



\*i) = Gf(E)Vi\* l ) 



(3.4) 



at real negative and complex E energies, respectively. 

We are going to solve these equations in a unified way 
by approximating only the potential term Vi. For this 
purpose we write the unit operator in the form 



li = lim l lN 

N— >oo 



where 
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The a factors have the properties \mi n ^^ L a n 



N 



(3.5) 



(3.6) 



and 



limjv— >oo o - " = 1, and make the limit in ( |3.5| ) smoother. 
They were introduced originally for improving the con- 
vergence properties of truncated trigonometric series (lB) , 



Vi = UVlh W UnViUn = V, N = ^ I") ^nn' 



where the matrix elements 

V nn ,=a%{n\Vi\n')a» , 



(3.8) 



(3.9) 



in general, have to be calculated numerically. This ap- 
proximation is called separable expansion, because the 
operator , e.g. in coordinate representation, appears 
in the form 



JV 



(r\V N \r')= <"V> - 



(3.10) 
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i.e. the dependence on r and r' appears in a separated 
functional form. 



With this separable potential Eqs. (3.3) and (3^4) re- 
duce to 



N 



*r> = |$f ) )+ £ G?(E±iO)\n)V nn ,(n'\V 



(±)\ 
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and 



|¥,> = £ G?{E)\n)V nn , (n'|¥, 



(3.11) 



(3.12) 
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respectively. To derive equations for the coefficients 
vj/ ; (±) = (n'l*^) and = (n'\^i), we have to act 
with states (n"\ from the left. Then the following in- 
homogeneous and homogeneous algebraic equations are 
obtained, for scattering and bound-state problems, re- 
spectively: 



[(G?(E±iO))-i-V l ]*< i ±) 



(±) 



and 



[(G?(E))- 1 -Vi]% = 0, 



(3.13) 



(3.14) 



where the overlap $ij = (n|$j ) can also be calcu- 



lated analytically [|4j|. The homogeneous equation (3.14) 
is solvable if and only if 
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det[(Gp(25))- 



V, 







(3.15) 



holds, which is an implicit nonlinear equation for the 
bound- and resonant-state energies. As far as the scat- 
tering states concerned the solution of ( |3.13 ) provides 
the overlap (n|\&;). From this quantity any scattering 
information can be inferred, for example the scattering 
amplitude corresponding to potential V\ is given by Jl^ ] 



A 



* I ( - ) i^i*ro=* 



i(-) 



(+) 



(3.16) 



Note that also the Green's matrix of the total Hamilto- 
nian, which is equivalent to the complete solution of the 
phys ical system, can be constructed as the solution of 
Eq. 



G l (z) = [(G?(z)y 



■n 



(3.17) 



Finally, it should also be emphasized, that in this ap- 
proach only the potential operator is approximated, but 
the asymptotically important H c term remains intact. 
The solutions are defined on the whole Hilbert space, 
not only on a finite subspace. The wave functions are 
not li near c ombin ation s of basis functions, but rather, as 
Eqs. (3.11) and (3.12) indicate, linear combinations of 



the states Gf(E)\n), which have been shown to possess 
correct Coulomb asymptotics [||. 



A. Bound, resonant and scattering states in a model 
nuclear potential 

In this subsection we apply the above techniques to 
calculate bound-, resonant- and scattering-state solutions 
of potential problems. The particular example we con- 
sider is a potential representing the interaction of two 
a particles. This example is thoroughly discussed in 
the pedagogical work |l7J in the context of a conven- 
tional approach based on the numerical solution of the 
Schrodinger equation. The interaction of two a particles 
can be approximated by the potential 



V a - a (r) = -Aexp(-f3r 2 ) + 



Z 2 c 2 



erf( 7 r), (3.18) 



where erf(z) is the error function [13]. This potential 
is a composition of a bell-shaped deep, attractive nu- 
clear potential, and a repulsive electrostatic field between 
two extended charged objects. The units used in the 
Hamiltonian of this system are suited to nuclear phys- 
ical applications, i.e. the energy and length scale are 
measured in MeV and fm, respectively. In these units 
h/(2m) — 10.375 MeV fm 2 (with m being the reduced 
mass of two a particles) and e 2 = 1.44 MeV fm. The 
other parameters are A =122.694 MeV, = 0.22 fm" 2 , 
7 = 0.75 fm -1 and Z = 2 (the charge number of the a 
particles) . 

Since the a — a potential possesses a Coulomb tail, 
for the asymptotically relevant operator we should take 



the Coulomb Hamiltonian H c of Eq. ( |2.l[ ) and the short- 
range potential is then defined by 

Z 2 e 2 Z 2 e 2 

V(r) = V a - a (r) = -Aexp(-/?r 2 ) erfc(7r), 

r r 

(3.19) 

with erfc(z) = 1 — erf(z). We approxi mat e the short- 
range potential V(r) according to Eq. ( |3.8| ) on the CS 
basis. From this point on, however, the properties of the 
physical system will be buried into the numerical values 
of the matrix elements. Thus the method is applicable to 
all types of potentials, as long as we can calculate their 
matrix elements somehow. Besides ordinary potentials 
this equally applies to complex, momentum-dependent, 
non-local, etc. potentials relevant to practical problems 
of atomic, nuclear and particle physics. Furthermore, the 
present formalism is equally suited to problems includ- 
ing attractive or repulsive long-range Coulomb-like and 
short-range potentials. 

It should also be emphasized again that in this ap- 
proach only the potential operator V is approximated 
and the resulting wave function possesses the correct 
asymptotics. So, if the parameter b is chosen such that 
V in (3.£) approaches V uniformly, the convergence 
of the bound-, resonant- and scattering-state quantities 
with respect to increasing TV will also be uniform, thus 
V N provides more or less uniformly good approximation 
to V on the whole spectrum of physical interest. 



1. Bound states 

First we consider only the nuclear part of potential 
( p. 18 ) and switch off the Coulomb interaction by setting 
Z = 0. According to Rcf. |17|] , this potential supports 
altogether four bound states: three with I — and one 
with I = 2. However, it is known that the first two I = 
and the single I = 2 state are unphysical, because they 
are forbidden due to the Pauli principle. This fact is 
not taken into account in this simple potential model. 
Although from the physical point of view these Pauli- 
forbidden states have to be dismissed as unphysical, they 
are legitimate solutions of our simple model potential. 
The proper inclusion of the Pauli principle in the model 
would turn the potential into a non-local one. This prob- 
lem has been considered within the present method in 
Ref. §. 

As illustrative examples we present the results of our 
calculations for the three I = stat es. W e determined 
the energies of these states from Eq. ( |3.15 ), using the CS 
parameter b = 4 fm -1 . Table [| shows the convergence 
of the method with respect to N, the number of basis 
states used in the expansion. It can be seen that the 
method is very accurate, convergence up to 12 digits can 
easily been reached. We note that according to Ref. |jl7) , 
the energy of the two lowest (i.e. the unphysical) I = 
states is E = -76.903 6145 and -29.000 48 MeV in the 
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uncharged case, which is in reasonable agreement with 
our results. 



2. Resonance states 

Switching on the repulsive Coulomb interaction (Z = 
2) the bound states are shifted to higher energies. The 
most spectacular effect is that the third I = state, which 
is located at E = -1.608 740 8214 MeV in the uncharged 
case, moves to positive energies and becomes a resonant 
state. This is in agreement with the observations: the 
a — a system (i.e. the 8 Be nucleus) does not have a stable 
ground state, rather it decays with a half life of 7 x 10~ 17 
sec. 

In our calculations we determined the energy corre- 
sponding to this resonance and to other ones as well by 
the same techniques we used before to find bound states. 
In fact, we used the same computer code and the same 
CS parameter (b — 4 fm _1 ) as we used in the analysis of 
bound states. The method, again, requires locating the 
poles of the Green's matrix, but not on the real energy 
axis, rather on the complex energy plane. In Table [n] we 
demonstrate the convergence of our method with respect 
to N for the lowest I = and I = 2 resonance states. In 
Fig. H we plotted the mo dulus of the determinant of the 
Green's matrix (as in Eq. (3.17)) \G(E)\ over the complex 
energy plane for 1 = 2. The resonance is located at the 
pole of this function. Finally, we mention that there is a 
resonance state for I = 4 at E = 11.791 038 - i 1.788 957 
MeV. 

Although it is not our aim here to reproduce ex- 
perimental data with this simple potential model, we 
note that the corresponding experimental values ]isf| 
are E Td=0 =0.09189 MeV, E^ i=2 =3.132 ± 0.030 MeV, 
# r ,i=4 =11.5 ± 0.3 MeV, and r /=0 /2 = (3.4 ±0.9) x 10~ 6 
MeV, r ;=2 /2=0.750 ± 0.010 MeV, r /=4 /2 -1.75 MeV. 



3. Scattering states 

In order to demonstrate the performance of our ap- 
proach for scattering states we calcu lated scattering 
phase shifts b"i{E) for V a - a (r) in ( 3.18j ). As described 
previously in this section, phase shifts can be extracted 
from the scattering amplitude given in Eq. (3.16). Spec- 



ifying this formula for the Coulomb-like case and for a 
given partial wave I we have 



ai = - exp(i(2?7 ; + Si)) sin Si , 



(3.20) 



where a/ is the Coulomb-modified scattering amplitude 
corresponding to the short-range potential, rji = argr(Z + 
ij+l)) is the phase shift of the Coulomb scattering with 
7 = Z 2 e 2 m/h 2 k and Si is the phase shift due to the 
short-range potential. 



The convergence of the method with respect to N is 
demonstrated in Table III, where Sq(E) is displayed at 
three different energy values E. As in our calculations 
for the bound and the resonance states, we used 6 = 4 
fm _1 here too. 

In Fig. H we plotted the scattering phase shifts Si (E) for 
I = 0, 2 and 4 up to E = 30 MeV. In all three plots in Fig. 
^ the location of the corresponding resonance is clearly 
visible as a sharp rise of the phase around the resonance 
energy E T . This rise is expected to be more sudden for 
sharp resonances, and this is, in fact, the case here too. 
The phase changes with an abrupt jump of tt for the 
sharp I = resonance, while it is slower for the broader 
I = 2 and I = 4 resonances. We also note that the phase 
shifts plotted in Fig. || are also in accordance with the 
Levinson theorem, which states that <5/(0) = to7t, where 
m is the number of bound states in the particular angular 
momentum channel. Indeed, as we have discussed earlier, 
there are two bound states for I = 0, one for I = 2 and 
none for I = 4. 

As an illustration of the importance of the smoothing 
factors we show in Fig. [I] the convergence of the phase 
shift Sq(E) at a specific energy E = 10 MeV with and 
without the sm oothing factors in (3.6), and conse- 
quently in (|3.9|). (Here and everywhere else the a pa- 



rameter was chosen to be 5.2.) Clearly, the convergence 
is much poorer without the smoothing factors. We note 
that this also applies to the other quantities calculated 
for bound and resonance states. 

Finally, we should call the attention upon the fact that 
this method is very accurate. Reasonable accuracy is 
reached already at relatively small basis, around N = 20. 
The accuracy gained in larger bases is beyond most of the 
practical requirements. Test calculations have been per- 
formed on a linux PC (Intel PII, 266 MHz) using double 
precision arithmetic. The calculation of a typical bound- 
or resonant-state energy requires the evaluation of the po- 
tential matrix by Gauss-Lague rre qu adrature and finding 
the zeros of the determinant ( 3.15 ), which incorporates 
the evaluation of the Coulomb Green's matrix and the 
the calculation of a determinant by performing an LU 
decomposition in each steps. The determination of the 
energy value in the first column and last row of Table 
H which meant 6 steps in the zero search and handling 
of 40 x 40 matrices, took 0.06 sec. The corresponding 
resonance energy value in Table || required 12 steps in 
the zero search on the complex energy plane and 0.8 sec. 
The evaluation of the three phase shift values in the last 
row of Table III took 0.19 sec, 0.11 sec. and 0.08 sec, 



respectively. So, this method is not only extremely accu- 
rate but also very fast. 



IV. CONCLUSIONS 

In this work we have presented a continued frac- 
tion representation of the Coulomb Green's operator on 
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Coulomb-Sturmian basis. Numerical illustrations show 
that this representations is simple, readily computable 
and numerically exact on the whole complex energy 
plane. This result is based on the observation that in 
this basis the Coulomb Hamiltonian possesses a Jacobi 
matrix structure. These techniques can be transferred to 
other problems, where the Hamiltonian matrix also has 
a Jacobi form. Examples for this are the harmonic os- 
cillator ||, the generalized Coulomb problem Q, which 
contains both the D dimensional harmonic oscillator and 
the Coulomb problems as special limits and the relativis- 
tic Coulomb problem 

Once the representation of the Coulomb Green's oper- 
ator in the discrete CS basis is available, we can proceed 
to solve the Lippmann-Schwinger integral equation. In 
practice this means the approximation of the potential 
term on a finite subset of this basis. This is the only stage 
where approximations are made, otherwise this method is 
exact and analytic, and provides asymptotically correct 
solutions. Consequently, bound, resonance and scatter- 
ing problems can be treated on an equal footing, while 
these phenomena are usually discussed in rather different 
ways in conventional quantum mechanical approaches. 
This unified treatment is also reflected by the fact that 
all the calculations are made using the same discrete ba- 
sis, containing also the same basis and other parameters. 
Furthermore, this approach is applicable to a wide range 
of potential problems, including complex, momentum- 
dependent, non-local, etc. ones. These techniques were 
illustrated with calculations for a model nuclear potential 
describing the interaction of two a particles. 

Once a discrete Hilbert-space basis representation of 
the Green's operator of the free or Coulomb Hamilto- 
nian is at our disposal, we can also construct the ma- 
trix representation of the Green's operator of the to- 
tal Hamiltonian by solving a resolvent equation. This 
provides a complete description of the two-body system. 
The Green's operators, similarly to the wave functions, 
contain all the information about the system, and any 
physical property can be derived from them. The impor- 
tance of this result is not purely formal, but it provides us 
with an effective tool in practical calculations for realistic 
physical problems. Moreover, from the Green's matrices 
of simpler systems, the Green's matrices of more complex 
composite systems can be constructed. These techniques 
have been used so far for solving Faddeev equations of 
some three-body Coulombic systems |t],|8| but they can 
certainly be adapted to other challenging problems as 
well. 
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TABLE I. Convergence of the I = bound-state energy 
eigenvalues E n i in V a - a (r) in the uncharged (Z = 0) case. N 
denotes the number of basis states used in the expansion. 



G 



N E 00 (MeV) 



Bio (MeV) 



E 2() (MeV) 



8 


-76.903 


557 


1529 


-29.005 


234 


9134 


-1.739 


478 


2626 




10 


-76.903 


609 


9717 


-29.000 


352 


3141 


-1.637 


269 


0831 




15 


-76.903 


614 


3090 


-29.000 


469 


8249 


-1.608 


824 


6403 




18 


-76.903 


614 


3254 


-29.000 


470 


2338 


-1.608 


742 


5166 




20 


-76.903 


614 


3263 


-29.000 


470 


2566 


-1.608 


741 


0685 


LLJ^ 


25 


-76.903 


614 


3265 


-29.000 


470 


2623 


-1.608 


740 


8256 


O 


28 


-76.903 


614 


3265 


-29.000 


470 


2625 


-1.608 


740 


8216 


O 


30 


-76.903 


614 


3265 


-29.000 


470 


2625 


-1.608 


740 


8213 


CD 


35 


-76.903 


614 


3265 


-29.000 


470 


2626 


-1.608 


740 


8214 




40 


-76.903 


614 


3265 


-29.000 


470 


2626 


-1.608 


740 


8214 





TABLE II. Convergence of the energy eigenvalues E TeBi i 
for the I = and I = 2 resonances in the V a - a (r) potential. 
N denotes the number of basis states used in the expansion. 



N 






3 re s,0 


(MeV) 






E 


res, 2 


(MeV) 




8 


-0.000 


854 


9596 +i 0.000 


000 


0000 


2.807 


21 


-i 0.607 


11 


10 


0.063 


364 


2503 


-i 0.000 


000 


0681 


2.866 


30 


-i 0.628 


56 


15 


0.091 


785 


0787 


-i 0.000 


002 


8092 


2.889 


68 


-i 0.620 


99 


18 


0.091 


963 


0277 


-i 0.000 


002 


8572 


2.889 


34 


-i 0.620 


53 


20 


0.091 


969 


7296 


-i 0.000 


002 


8588 


2.889 


24 


-i 0.620 


56 


25 


0.091 


971 


8479 


-i 0.000 


002 


8592 


2.889 


23 


-i 0.620 


62 


28 


0.091 


971 


9788 


-i 0.000 


002 


8592 


2.889 


25 


-i 0.620 


62 


30 


0.091 


972 


0064 


-i 0.000 


002 


8592 


2.889 


25 


-i 0.620 


61 


35 


0.091 


972 


0258 


-i 0.000 


002 


8592 


2.889 


24 


-i 0.620 


61 


40 


0.091 


972 


0290 


-i 0.000 


002 


8592 


2.889 


25 


-i 0.620 


61 



-0.2 







-0.15 -0.1 -0.05 
E (atomic units) 

FIG. 1. The determinant of a 3 x 3 Coulomb Green's ma- 
trix det\G c (E)] as the function of the energy E for / = 1. The 
bound states of the Coulomb problem are located at energies 
where the vertical lines cross the horizontal axis. (These lines 
are shown only for demonstrative purposes: they do not cor- 
respond to functional values of det[G c (£)].) Atomic units of 
m = e = h — 1 and Z = — 1 were used. For the sake of clarity 
only the first 6 energy levels are shown. These are located at 
E n = -l/[2(n + I + l)] 2 with n < 5. 



TABLE III. Convergence of the 80(E) phase shift (in ra- 
dians) in the V a - a {r) potential at three different energies. N 
denotes the number of basis states used in the expansion. 



N 


E = 0.1 MeV 


E = 1 MeV 


E = 30 MeV 


8 


6.283 230 


8.817 731 


7.783 217 


10 


9.424 059 


8.862 581 


4.817 163 


15 


9.424 018 


8.859 651 


4.835 479 


18 


9.424 022 


8.859 467 


4.829 861 


20 


9.424 023 


8.859 441 


4.828 882 


25 


9.424 024 


8.859 419 


4.828 563 


28 


9.424 024 


8.859 414 


4.828 555 


30 


9.424 024 


8.859 412 


4.828 554 


35 


9.424 024 


8.859 411 


4.828 552 


40 


9.424 024 


8.859 411 


4.828 552 



■ 

--Si 




FIG. 2. The modulus of the determinant of the G(E) 
Green's matrix for the a — a potential on the complex en- 
ergy plane for I = 2. The pole at E = 2.889 25 - i 0.620 61 
MeV corresponds to a resonance. 
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10 20 30 

E(MeV) 

FIG. 3. Scattering phase shifts Si(E) (in radians) in the 
a — a potential for I — 0, 2 and 4. The resonances in these 
partial waves appear as sharp rises in the corresponding phase 




N 

FIG. 4. Convergence of the scattering phase shift 5o(E) 
(in radians) at E — 10 MeV calculated with (solid line) and 
without (dashed line) smoothing factors. 
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